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SUMMARY 

A database for YO, ,-NdO^-YbO, 5 -Sc0 i5 -Zr0 2 for ThermoCalc* has been developed. The basis of this work is 
the YO J5 -ZrO, assessment by Du et al. (ref. 1). Experimentally only the YO, s -Zr0 2 system has been well-studied. 
All other systems are only approximately known. The major simplification in this work is the treatment of each 
single cation unit as a component. The pure liquid oxides are taken as reference states and two term lattice stability 
descriptions are used for each of the components. The limited experimental phase diagrams are reproduced. 


I. INTRODUCTION 

The rare earth oxide stabilized zirconias are of considerable technological importance. Applications range from 
thermal barrier coatings for turbine blades to electrodes for sensors and fuel cells. The most widely used system is 
Y 2 0 3 -Zr0 2 . The phase diagram for this system has been extensively studied (refs. 2 to 11) Further, there are several 
thermodynamic assessments and calculated phase diagrams in the literature for this and related systems (refs. 1,12 
to 19). 

In this report we also consider the stabilizers: Yb 2 0 3 , Sc 2 0 3 , and Nd 2 0 3 . These different stabilizers alone and in 
combinations may offer improved properties. However, both experimental phase diagrams and thermodynamic data 
for these systems are very limited (refs. 20 to 30). Reliable phase and thermodynamic information over a range of 
compositions and temperatures is essential to understand the processing and properties of zirconia w ith Y 2 0 3 , Yb 2 0 3 , 
Sc 2 0 3 , and Nd,0 3 stabilizers. 

Figures 1 to 10 give the available psuedo-binary phase diagrams for the systems (refs. 22 to 30) listed in 
table I. The most complete and most recent phase diagrams are given. In the case of Y 2 0 3 -Zr0 2 , two phase diagrams 
are shown in figures 1(a) and (b). The first diagram from Stubican et al. (refs. 1 1 and 22) includes one intermediate 
compound— Zr 3 Y 4 O l2 . The second diagram (ref. 23) also includes an ordered hexagonal phase— ZrY 6 0 (2 . However, 
the existence of the second compound remains controversial and will not be considered here. In addition, the pri- 
mary interest is in the Zr0 2 -rich portion of the diagram and therefore its existence is not important in the present 
study. Figures 2(a) and (b) are the available diagrams for the Yb,0 3 -ZrO, system. Figure 2(a) primarily shows the 
liquidus and figure 2(b) shows the solid phases. Note that these are similar to the Y 2 0 3 -Zr0 2 diagram and shows an 
analogous intermediate phase — Zr 3 Yb 4 O t2 . How ever Zr 3 Yb 4 0, 2 is stable to 1612 °C and Zr 3 Y 4 O l2 is stable to only to 
1250 °C (fig. 1(a)) or 1330 °C (fig. 1(b)). Figures 3 to 10 are the remaining available binaries. Most of these dia- 
grams are incomplete and have many approximate boundaries. ......... 

The addition of two (or more) stabilizing rare earth oxides to zirconia may lead to some interesting properties. 
Yet only a portion of the Y 2 0 3 -Yb,0 3 -Zr0 2 ternary diagram has been determined experimentally (ref. 31). Three 
isothermal sections are shown in figures 1 1(a) to (c). The major features of these systems are a continuous solid so- 
lution of tetragonal zirconia and cubic zirconia with both Y 2 0 3 and Yb 2 0 3 . Note also that the two intermediate com- 
pounds — Zr 3 Y 4 O l2 and Zr 3 Yb 4 O i2 connect at lower temperatures. At higher temperatures Zr 3 Y 4 0 12 becomes unstable 
and only Zr 3 Yb 4 O p protrudes into the ternary. At the highest temperatures, both compounds are unstable. No other 
ZrO, rare earth oxide experimental ternaries were found. 

The first calculations of these types of refractory oxide phase diagrams were done in 1988 by Kaufman 
(ref. 14). Kaufman treats each oxide unit as a component. He sets the liquid oxide as a reference state and uses two- 
term lattice stabilities to describe each end point oxide. Thus the lattice stability of a particular polymorph becomes 
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the energy to form that polymorph from the liquid. Du et al., (refs. 1 and 15) have done two assessments of the 
Y,0 3 -Zr0 2 , using the Lukas programs (ref. 32) to derive phase descriptions that best fit all the available experimental 
data. Their first assessment (ref. 1) also uses the liquid as the reference state. Their second assessment (ref. 15) in- 
vokes the standard element reference state (SER) and uses Gibbs energy expressions for each phase. 

This report discusses the development of a Thermo-Calc database for the multicomponent system: Nd,0 3 -Sc 2 0 3 - 
Yb 2 0 3 -Y,0 3 -Zr0 2 . The quantities of Du et al., for Zr0 2 -Y,0 3 (ref. 1) are used directly and other binaries are devel- 
oped based on this. Du et al., (ref. 1) perform two optimizations for the ZrO,-Y,0 3 system — using a different set of 
liquidus data for each optimization. Their second optimization values are used here. The rare earth oxides are nor- 
malized to one cation: Nd0 15 -Sc0 l 5 -Yb0, ,-YO L 5 -ZrO r Each of these one cation oxide units are treated as a compo- 
nent. Two-term lattice stability phase descriptions referenced to the liquid state are used, which are less complex 
than the Gibbs energy phase descriptions based on the standard element reference state. 

The entire database in standard Thermo-Calc *.tdb format is given in appendix I. In the following report, we 
discuss this database and the calculation of binary and higher order diagrams in the Nd 2 0 3 -Sc 2 0 3 -Yb,0 3 -Y 2 0 3 -Zr0, 
system. 


II, LATTICE STABILITIES FOR END-POINTS 

In order to estimate the lattice stabilities with respect to the liquid state, we need the enthalpies and entropies for 
the polymorphic transformations. Only limited information is available, particularly on the solid/solid transforma- 
tions. These data (refs. 33 to 35) are shown in table II. As noted, each oxide is normalized to one mole of cation. The 
oxide abbreviations are given in table III. 

It is important to discuss the compositional difference between the calculations and the measured values in 
table II and figures 1 to 1 1. The experimental phase diagrams are for RE,0 3 -ZrO 2 (RE = Rare Earth Element) and 
will be referred to as such; the calculated diagrams are for REO, 5 -ZrO, and will be referred to as such. In order to 
compare the two, we need to convert from mole fraction x(RE 2 0 3 ), in RE 2 0 3 -Zr0 2 to mole fraction x(RE0 15 ), in 
REO i 5 -ZrO,. Note that 


x(RE 2 0 3 in RE 2 O 3 -Zr0 2 ) = 


n(RE 2 03 in RE 2 0 3 — Zr0 2 ) 

n(RE 2 0 3 in RE 2 0 3 -Zr0 2 ) + n(Zr0 2 in RE 2 0 3 -Zr0 2 ) 


( 1 ) 


Here n is the number of moles. The composition for REO, 5 -Zr0 2 is: 

x(REO, 5 in REO! 5 - Zr0 2 ) = n( REO l.5 ln REQ !,5 ~ Zr ° 2 ) 

n(REOj 5 in REOJ 5 - Zr0 2 ) + n(Zr0 2 in REO, 5 - Zr 0 2 ) 


Note that for each composition n(Zr0 2 ) the same and that n (RE,0 3 ) = Vi n (REO, 5 ). So now: 


x(REO] 5 ) 


2x(RE 2 0 3 ) 

2x(RE 2 0 3 ) + (l -x(RE 2 0 3 )) 


( 2 ) 


( 3 ) 


Table IV shows the equivalent compositions. 

Now consider the derivation of the lattice stabilities for each oxide. The most information is available for Zr0 2 , 
which undergoes three transformations (ref. 36): 

Zr0 2 : Monochnic — » Tetragonal — > Cubic — > Liquid (4) 

The approach of Kaufman (ref. 14) and ThermoCalc notation are used in the following discussion. 

The liquid phase is taken as the reference state and thus: 

G(L, ZM; 0) = 0.0 (5) 
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For the cubic phase: 


G(C,ZM; 0) = AH C I - T AS CL = -87027 + 29.47 1 T (6) 

This is the value from Kubachewski and Alcock (ref. 33) with a melting point of 2953 K. More recent 
investigations have given the melting point of cubic zirconia as 2983 K. Du and Jin (ref. 1) adjust equation (6) to 
reflect this: 


G(C, ZM; 0) = -87986.6 + 29.496T (7) 

The tetragonal to cubic transformation temperature is given as 2568 K (33) and thus: 

AH tc -T AS t c = -5941 + 2.313T (8) 

Again Du and Jin (ref. 1) adjust this for a more recent transformation temperature of 2642 K: 

AH t c - T AS t c = -594 1 + 2.249T (9) 

In order to obtain the free energy of the tetragonal phase referenced to the liquid, equations (7) and (9) are 
added: 


G(T, ZM; 0) = AH tt - T AS TI = -93927 + 3 1 .745T ( 1 0) 

From equation (10) a ‘melting point’ for the tetragonal phase can be calculated to be 2959 K. A metastable 
melting point must always be less than the melting point of the stable phase. 

Table II indicates some disagreement on the monoclinic to tetragonal transformation temperature and heat. Here 
we take the temperature as 1454 K and the heat as 6000 joules/mol: 

AH mt - T AS M T = -6000 + 4. 1 27T (11) 

The free energy of the monoclinic phase referenced to the liquid is obtained by adding equations (10) and (11): 

G(M, ZM; 0) = AH mi - T AS ML = -99927 + 35.872T (12) 

For the other oxides — Y0 15 , YbO, 5 , Nd0 15 , ScO, -a similar approach was followed. However, less data is 
available. Of these the most data is available for YO i5 . Y0 15 undergoes two phase transformations: 

Y0 15 : Cubic Hexagonal — > Liquid (13) 

Again, the liquid phase is taken as the reference state: 

G(L, YM; 0) = 0.0 (14) 

The melting point of the hexagonal phase is accepted as 2712 K. Kaufman (ref. 14) estimates the lattice stability 
of the hexagonal phase as: 


G(H, YM; 0) = AH h L - T AS HL = -56735 + 20.92T ( 1 5) 

This enthaply is somewhat more than half the IVTAN (ref. 34) data for melting of Y 2 0 3 , but more in line with 
Kubachewski’s (ref. 34) values for alumina heat of melting (54 kJ/mol AlO, 5 ). 

The only measured transition heat for the polymorphs of a rare earth oxide was that for Y 2 0 3 cubic to hexago- 
nal, as reported in IVTAN to be 54 kJ/mol. However, this number seems much too high for a solid/solid transition 
such as this. Kaufman (ref. 14) estimates the heat to be about 10 kJ/mol and it seems more reasonable to use this 
value. Accepting the cubic to hexagonal transformation temperature of 2550 K the lattice stability for the cubic 
phase becomes: 


G(C,YM; 0) = AH c l - T AS CL = -67419 + 25.105T 


(16) 
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The behavior of YbO, 5 is controversial. The two published phase diagrams with YbO i3 do not agree. The 
Yb0 i 5 -Zr0 2 diagram (fig. 2) indicates that the hexagonal phase melts; whereas the YbO l5 -NdO i5 diagram (fig. 9) 
indicates that the cubic phase melts. It is assumed that the hexagonal phase of YbO J5 melts, in analog} 7 with YO i5 . 
Therefore the phase sequence with temperature is taken as: 

YbOj 5 : Cubic — » Hexagonal — > Liquid 

For NdO, the phase sequence with temperature appears to be: 

NdO, 5 : N phase — » Hexagonal — > Cubic — * Liquid 

Little information is available for ScO, . and it is assumed that the cubic phase melts. From the melting points 
and estimated heats of fusion, the lattice stabilities are derived as shown in the database (appendix I) 


III. INTERACTION PARAMETERS: CUBIC, MONOCLINIC, TETRAGONAL, N PHASE, NM PHASE 


In the binary case, each of the above solution phases are described by the following expression for the free 
energy: 


2 2 

G m (x,T) = ^x i 0 G i (D+RT^Xj lnx; + ex G m (17) 

i=l i=l 

The first term is the sum of the lattice stabilities for each of two components. These are discussed in the previ- 
ous section. The Calphad approach requires a lattice stability for every component in each phase — even if a particu- 
lar component does not form the phase described. For example, a lattice stability if necessary for hexagonal Zr0 2 , 
monoclinic YO ]5 , and tetragonal YO } 5 . The lattice stabilities are chosen so that these phases are always unstable 
(Appendix I). 

The second term in equation (17; is the random mixing term for an ideal solution. The third term is the excess 
Gibbs energy and accounts for the deviations from ideality. Many solution models have been developed for this term 
(refs. 37 to 39). In the Calphad approach, the excess Gibbs energy is expanded as a Redlich-Kister polynomial 
(ref. 37): 


ex G m =x 1 x 2 ^ Jl (x 1 -x 2 ) j 
j=0 


(18) 


Here only the first two terms are considered: 

ex G m = x x x^ L + x x x, L(x v -x 2 ) = x l x 1 °L + x i x 2 1 L(l-2x 2 ) ( 19 ) 


where °L and 'L are the zero and first order interaction parameters, respectively. In assigning mole fraction designa- 
tions, ThermoCalc takes the components in alphabetical order. Thus in this case x 2 = x ZK) and the sign of the second 
term in equation (19) is always the opposite of the sign of X for x 2 > 0.5. 

The shape of the Gibbs energy curves versus composition and temperature determine the stable phase fields in a 
classical temperature-composition phase diagram. The first two terms in equation (17) are fixed and lead to an ideal 
solution, as discussed in many textbooks (refs. 38 to 39). The variables in the excess term (X and X) are adjusted to 
fit experimental observations for a real solution, which generally show deviations from ideality. 

The YO, 5 -Zr0 2 system description can be simplified by introducing a miscibility gap in the cubic phase field. 
This will be discussed in the next section. It follow s that theX and X parameters need to be adjusted to create this 
feature. This is illustrated in figure 12. A single zero order term (regular solution) can only yield a symmetric misci- 
bility gap. The X term is also needed to produce an asymmetric miscibility gap. The effect of the sign of X on the 
location of the miscibility gap is shown in figure 12 and table V. For the cubic phase in a RE0 t 5 -Zr0 2 phase dia- 
gram, we want the miscibility gap on the Zr0 2 -rich side. This requires a positive X and a negative second term in 
equation (19), which gives a second minima on the right, but less than the first minima on the left, as shown in 
figure 12. 
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The interaction parameters for the five solid solutions in this system are thus determined by (a) Using the 
assessment of Du (1) as a guide and (b) Estimating numbers which give a G m which will lead to phase regions simi- 
lar to those observed in the experimental phase diagram. 


III. A. Cubic Solid Solution (Cubic ZrO,ss and Cubic RE0 13 ss) 

A major feature of the Zr0 2 -Rare Earth Oxide phase diagrams is the large cubic phase field. This is evident in 
each of the experimental phase diagrams (figs. 1 to 4). 

As pointed out by Degtyarev and Voronin (refs. 1 2 and 1 3), most of the Y,0 3 -Zr0 2 phase diagram is composed 
of cubic solid solutions of Y0 15 and ZrO,. Cubic zirconia has a classical face-centered cubic lattice CaF. structure. 
The cubic phase for Y,0 ( has a body-centered cubic lattice like Mn 2 0 3 . Degtyarev and Voronin (refs. 12 and 13) 
further point out that this Mn,0 3 structure can be derived from the CaF 2 structure by removing one quarter of the 
oxygen anions. Thus the two cubic phases can be derived from one G-x curve with a miscibility gap. 

The above assumption simplifies the phase diagram calculation procedure. The miscibility gap option in 
Thermo-Calc is used which has the program search for a second free energy minimum. Further, the experimental 
phase diagrams indicate that the miscibility gap is shifted to the right. 

The interaction parameters for the YO, s -ZrO, are taken from the assessment of Du et al. (ref. 1): 

L(C,YM,ZM; 0) = -12060 + 1 1. 156T (20) 

L(C,YM,ZM; I) = +13784 + 5.379T (21) 

The first order interaction parameter is always positive and thus shifts the miscibility gap to the right (table V) 
as observed experimentally. 

As shown in the database (appendix I), similar interaction parameters were taken for the YbO, 5 -Zr0 2 and ScO l5 - 
ZrO, systems. For the NdO, 5 -ZrO, system, the following interaction parameters are estimated: 

L(C, NM, ZM; 0) = -52400 + 27T (22) 

L(C,NM,ZM; 1) = +24400 +1.5T (23) 

Again the first order interaction parameter is always positive and hence the miscibility gap is shifted to the right, as 
the experimental phase diagram shows. Note that equation (23) yields the similar values to equation (21) at high 
temperatures. 

Consider the remaining interaction parameters between the various rare earth oxides. Now there is no miscibil- 
ity gap and in most instances a continuous or nearly continuous cubic solution (figs. 5 to 10). Only the Yb 2 0 3 -Nd 2 0 3 
and Y,0 3 -Nd 2 0 3 systems have been established with any certainty. These interaction parameters are estimated as: 


L(C, BM,NM;0) = +21425 + 2.5T (24) 

L(C,BM, NM; 1)= -19600 + 8T (25) 

L(C,NM, YM;0) = +21425 + 2.5T (26) 

L(C, NM, YM; 1) = -19600 + 8T (27) 


The other rare earth oxide phase diagrams (Sc 2 0 3 -Y,0 3 , Sc 2 0 3 -Nd 2 0 31 , Yb 2 0 3 -Sc 2 0 3 , Yb 2 0 3 -Y 2 0,) are only very 
approximately known (figs. 5, 6, 8, and 10, respectively). In this case only the zero order interaction parameter is 
estimated to be +20,000 J/mol. 


III.B. Hexagonal Solid Solution (Hex REO l5 ss) 

As shown in figures 1(a) and (b), the hexagonal solid solution only occupies a small region of the Y 2 0 3 -Zr0 2 
phase diagram. The assessment of Du et al., gives a single Redlich-Kister term: 
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L(H, YM, ZM; 0) = 50420 (28) 

As shown in figure 12, a large positive value of °L gives a symmetric free energy of mixing with two small min- 
ima on either side. The minima on the YO i5 side gives the small hexagonal region in the phase diagram. 

Most of the other hexagonal °L solution terms were estimated to 20000 or 35000 J/mol. However the phase dia- 
grams for the Nd 2 0 3 -Y a 0 3 and Nd,0 3 -Yb,0 3 have been reported (figs. 7 and 9, respectively). In order to obtain the 
hexagonal solid solutions across these phase diagrams, the following expressions are used: 


L(H, BM, NM; 0) = 12750 +1.5T (29) 

L(H, BM, NM; 1) = -1000 - T (30) 

L(H, YM, NM; 0) =11265 +1.5T (31) 

L(H, YM, NM; 1) = 1000 - T (32) 


Both these sets of values for the interaction parameters give a free energy with a single, broad minima at the 
temperatures of interest. 


III.C. Liquid Solution 

The Liquid solution parameters for the YO tr ZrO, system have been determined by Du et aL, (ref. 1) to be: 


L(LIQ, YM, ZM; 0) = -183751 + 72.4T (33) 

L(LIQ, YM, ZM; 1 ) = 48733 - 9.476T (34) 

Due to similarity between the liquid regions in the Y0 l5 -Zr0 2 and Yb0 t5 -Zr0 2 systems, these same parameters 
are taken: 

L(LIQ, BM, ZM; 0) = -183750 + 72.4T (35) 

L(LIQ, BM, ZM; 1 ) = 48700 - 9.48T (36) 

For the NdO l5 -ZrO, system, a single term expression is used: 

L(LIQ, NM, ZM; 0) = -20000 (37) 

L(LIQ, NM, ZM; 1) = -15000 (38) 


All of the ScO l5 parameters were estimated to 20000 J/mol. As shown in figures 3, 6, 8, and 10, virtually noth- 
ing is known about the liquid phase in the scandia systems. Thus, it is reasonable to use the above approximations. 
For the other systems, zero order parameters are estimated as: 

L(LIQ, BM, NM; 0) = -200 (39) 

LfLIQ, BM, YM; 0) = -10000 (40) 

L(LIQ, NM, YM; 0) = -200 (41) 


III.D. Monoclinic Phase Solid Solution (Mon Zr0 2 ss) 

The monoclinic phases show 7 very limited solubility for a second oxide, except for the Nd0 15 -YbO 15 and 
Nd0 15 -Y0 15 systems. Again, beginning with the data from the assessment of Du et al., (ref. 1): 

L(M, YM, ZM; 0) = -58223 + 98.126T (42) 
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The other interaction parameters are estimated. Only two terms are used for the NdO, r YbO LS and 
Nd0 15 -Y0, , systems: 

L(M, BM, NM; 0) = -5000 (43) 

L(M, BM, NM; 1) = 100 (44) 

L(M, NM, YM; 0) = -50000 (45) 

L(M,NM, YM; 1)= 1000 (46) 

The remaining zero order interaction parameters are all estimated to be 20000 J/mol. 

im.E. Phase MN Solid Solution (MN ss) 

Again, since only the YbO, 5 -NdO L5 and YO u -NdO, 5 phase diagrams are known with any certainty, these inter- 
action parameters are taken as: 

L(MN, BM, NM;0) = -160000 + 48T (47) 

L(MN, BM, NM; 1) = 120000 - 15T (48) 

L(MN, NM, YM; 0) = -160000 + 50T (49) 

L(MN, NM, YM: 1) = 120000 + 15T (50) 

All other interaction parameters are taken as 20000 J/mol. 

III.F. N Phase Solid Solution (N ss) 

As above, since only the Yb0 15 -Nd0 l 5 and YO u -NdO, 5 phase diagrams are known with any certainty, these 
interaction parameters are taken as: 

L(N, BM, NM; 0) = -15200 -7T (51) 

L(N, BM, NM; 1) = -13000 - 6.7T (52) 

L(N,NM, YM;0) = -15000 -7T (53) 

L(N, NM, YM; 1) = 13000 + 5T (54) 

All other interaction parameters are taken as 20000 J/mol, except for: 

L(N, NM, ZM; 0) = 50000 (55) 

III.G. Tetragonal Phase Solid Solution (Tet Zr0 2 ss) 

This phase shows a triangular appearance in the Y 2 0 3 -Zr0, phase diagram (figs. 1(a) and (b)). There is a sug- 
gestion on the other RE 2 0 3 -ZrO, phase diagrams (figs. 2 to 4) of similar behavior. The assessment of Du et al., (1) 


L(T, YM, ZM; 0) = -25800 
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The same value is adopted for Yb0 15 -Zr0 15 : 

L(T, BM, ZM; 0) = -25800 (57) 

All other interaction parameters are estimated to be 20000 J/mol. 


IV. INTERMEDIATE COMPOUNDS 

In the binary NdO t5 -ZrO,, ScO (J -Zr0 2 , and Yb0 i5 -Zr0 2 systems there are several intermediate compounds. The 
Zr 3 Y 4 0 !2 phase has been reported by several investigators. These line compounds are treated with two sublattice 
model — the ZrO, component is on one sublattice and the Y0 15 component is one the other sublattice. Following 
Du et al.: 


T A ' 

G (J/7 )z 5 Y 4 Op = - °G C-Zr ° 2 +- °G C ~ Zr °2 - 3284.3- 2.26749T= -79536+ 24.72T (58) 

As mentioned, there is an analogous Zr 3 Yb 4 O l2 phase. This is treated as: 

G (I/7)Zr 3 Yb 4 0i2 = -79350 -1- 25T (59) 

To account for interchange of Y and Yb in this compound, the following solution parameter is introduced: 

L(DELTA,BM,YM:ZM;0) = 10200 (60) 

The free energy of each of the compounds are chosen so that they protrude into the cubic phase. For 
NdO, r Zr0 2 , there exists a binary with the 1:1 composition 

G(P, NM: ZM; 0) = -180450 + 57T (61) 

For the ScO, 5 -Zr0 2 , there are several intermediate compounds, each with a range of solution according to the 
experimental diagram in figure 3. However for the purpose of this approximation, we shall treat them as line com- 
pounds. 

For 4(Sc0 15 ) *3(Zr0 2 ), Phase RB: 

G(RB,SM:ZM;0)= -557200+ 171.5T (62) 

For 2(Sc0 15 ) ‘5(Zr0 2 ), Phase RH: 

G(RH, SM: ZM; 0) = -620000 + 1 85T (63) 

For 2(ScO I 5 ) *7(Zr0 2 ), Phase RM: 

G(RM, SM: ZM; 0) = -822800 + 258T (64) 


V. RESULTS: BINARY PHASE DIAGRAMS 

Having described the database, we now turn to the results for the binaries and compare them to the limited 
experimental data. Figures 13 to 22 give the calculated binaries. Recall again that the experimental diagrams are for 
RE 2 0 3 -Zr0, whereas the calculated diagrams are for REO, ZrO,. 

The calculated diagrams show the same general features as the experimental diagrams. The calculated 
Y0 J5 -Zr0 2 phase diagram (fig. 13) is, of course, the same as that obtained by Du et al., (ref. 1) and is in good 
agreement with the experimental diagram of Stubican et al., (fig. 1(a)). 

All other RE0 l s -ZrO, diagrams show 7 similar behavior to the Y0, 5 -Zr0 2 phase diagram on the Zr0 2 -rich side, 
where most of the interest Is. The NdO, 5 -Zr0 2 , ScO u -ZrO,, and Yb0 l5 -Zr0 2 diagrams are not so definitive, as much 
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less experimental data is available. Nonetheless, the general features of the experimental diagrams are reproduced on 
the calculated diagrams. 

The Yb 2 0 3 -Zr0, diagram is also only approximately known (fig. 2) and consists of a large cubic phase field and 
one intermediate compound. This is reproduced in the calculated diagram (fig. 14). 

The Nd,0 3 -Zr0 2 experimental diagram shows a eutectic at -21 15 °C and the NdO, 5 -ZrO, calculated diagram 
shows the same eutectic at 2100 °C. The Nd 2 0 3 -Zr0 2 experimental diagram shows a euctectoid at 1440 °C, whereas 
the calculated Nd0 15 -Zr0 2 phase diagram shows this eutectoid at -1660 °C. But this part of the diagram is not well- 
known. 

Only a limited portion of the Sc 2 0 3 -ZrO, phase diagram is known (fig. 3). For the purpose of this approxima- 
tion, the intermediate compounds (Sc 2 Zr 7 O l7 , Sc 2 Zr 5 O l3 , Sc.Zr.O,,) are treated as line compounds rather than solid 
solutions. Even with this approximation, the general appearance is reproduced. 

For the binaries among the rare earth oxides (Nd 2 0 3 , Sm,0 3 , Yb 2 0 3 , Y,0 3 ), only the Nd,0 3 -Y,0 3 (fig. 7) and 
Nd,0 3 -Yb,0 3 (fig. 9) systems are known to any extent. The major features of these diagrams are reproduced in the 
calculated diagrams (figs. 19 and 21, respectively). Both show a substantial monoclinic solid solution. The remain- 
ing rare earth oxide binaries show large cubic solid solution regions. 

There are some limited activity measurements for the Y 2 0 3 -Zr0 2 system at 2773 K (ref. 40) and the Sc 2 0 3 -Zr0 2 
system at 2600 K (ref. 41) obtained via mass spectrometry. In order to compare these data to calculations from the 
database described here, the mole fraction must be converted according to table IV and the activity also must be 
converted to appropriate units. Belov et al. (refs. 40 and 41) show that the rare earth oxides vaporize as: 

RE 2 0 3 (s) = 2 REO(g) + O(g) 

Thus the activity of RE,0 3 [a(RE 2 0 3 )] is defined as: 

a(RE 2 0 3 ) = ^ REO >ri P <-°JL 
[P°(REO)NP 0 (O)] 


(65) 


( 66 ) 


Here the quantities in the numerator are the partial pressures over the oxide in solution and the quantities in the 
denominator are the partial pressures over the pure oxide — solid Y,(), (presumably hexagonal) and cubic Sc.O,. 

Next consider RE0 1S , which should vaporize as: 

REO, 5 (s) = REO(g) + ‘/2 0(g) (67) 


a(REOj 5 ) 


[P°(RE0)][P°(0)] I/ - 


( 68 ) 


Thus the experimental measurements need only be converted with a square root for comparison to the calcu- 
lated values. 

Figures 23(a and b) and 24(a and b) show the comparison of these measured activities (refs. 40 and 41) to the 
calculated values from this database. For both systems the data indicate measurements were taken in a single phase 
region across the diagram; whereas the phase diagrams — both experiment and calculated indicate some two phase 
regions. There appears to be a conflict in the phase boundary locations. Agreement with a(Zr0 2 ) is good, as these 
data were used in the Du et al. (ref. 1) assessment. However, agreement with a(RE0 13 ) is only fair — these data were 
not used in the assessment. 


VI. RESULTS: TERNARY SECTIONS 

This database can be used to approximate ternary sections. One would expect the following features in a ternary 
diagram: 

1 . The cubic phase occupies a large region in each binary. Therefore it seems reasonable that the cubic phase 
would stretch across the ternary as well. 

2. The intermediate compounds should protrude into the ternary. 
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Calculations of the isothermal cuts led to some convergence problems in ThermoCalc and additional steps were 
necessary to force convergence. These are illustrated in the macro in appendix II. This macro also allows the varia- 
tion of the ternary interaction parameter. 

Figures 25(a) and (b) are isothermal cuts for the YO, 5 -YbO l5 -ZrO, system at 1600 K, illustrating the effects of 
two different ternary interaction parameters; 

Figure 24(a): L(C, BM, YM, ZM; 0) = +10000, insufficient for a continuous two phase region. 

Figure 24(b): L(C, BM, YM, ZM; 0) = +100000, continuous two phase region. 

In order to address issue (2) above, an interaction parameter was added for the intermediate phase. Figure 25 
shows the effect of 


L(C, BM, YM: ZM) = 70000 
L(DELTA, BM, YM: ZM) = 10200 

Figures 26(a) to (c) are the calculated ternary for 1473, 1673, and 2000 K. These can be compared to experi- 
mental sections in figure 1 1 . The protrusion behavior as a function of temperature for the delta phase is reproduced 
by these calculations. Further work is necessary to reproduce the phase boundaries of the solution phases correctly. 


VI. RESULTS: HIGHER ORDER SYSTEMS 

There are no reported data on higher order phase diagrams. However they can be approximated with this data- 
base. Phase data for systems with more than three components can be presented with a phase fraction diagram or an 
isoplethal section with one component varying over a small range. An example of a phase fraction diagram for 
ScO l5 -YO, 5 -YbO 15 -Zr0 2 is shown in figure 27. Figure 28 shows an isoplethal section for ScO i5 -YO i5 -YbO, 5 -ZrO, 
with a small variation in YO, 5 . 


VII. SUMMARY AND CONCLUSIONS 

An approximate database for ScO, 5 -Nd0 15 -Yb0 l5 -Y0 ]5 -Zr0 2 has been discussed. It is based on the assessment 
of Du et ah, for the Y0 t5 -Zr0 2 system (ref. 1) and uses a two-term lattice stability description of the phases with the 
liquid as the reference state. Binary phase diagrams for the other systems agree reasonably well with the limited 
experimental data. Ternary and higher order diagrams require the introduction of a ternary interaction parameter and 
an interaction parameter for the intermediate compounds to allow protrusion of the intermediate compounds into the 
multi-component composition space. The binary systems beyond Y 2 0 3 -ZrO, require additional experimental and 
modeling efforts. This will provide a stronger basis for further development of the multicomponent diagrams. 
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APPENDIX I 


$ Database file written 99-12-23 
$ 


ELEMENT 

/- 

ELECTRON_GAS 

0 . 0 00 0E + 0 0 

0 . 0000E+00 

0 . 0000E+00 

ELEMENT 

VA 

VACUUM 

0 . 0000E+00 

0 . 0000E+00 

0 . 0000E+00 

ELEMENT 

BM 

LI Q 

1 . 9704E+02 

0 . OOOOE+OO 

0 . 0000E+00 

ELEMENT 

NM 

LIQ 

1 . 6827E+02 

0 . OOOOE+OO 

0 . 0000E+00 

ELEMENT 

SM 

LIQ 

6 . 8960E+01 

0 . 0000E+00 

0 . 0000E+00 

ELEMENT 

YM 

LIQ 

1 * 1291E+02 

0 . 0000E+00 

0 . 0000E+00 

ELEMENT 

ZM 

LIQ 

1 .2322E+02 

0 . 0000E+00 

0 . 0000E+00 


FUNCTION UN_ASS 2.98140E+02 0.0 ; 3.00000E+02 N ! 

TYPE_DEFINITION % SEQ * ! 

DEFINE_SYSTEM_DEFAULT SPECIE 2 I 
DEFAULT COMMAND DEF SYS ELEMENT VA ! 


PHASE C % 1 1.0 ! 

CONSTITUENT C : BM, NM, SM, YM, ZM : ! 

PARAMETER G(C,BM;0) 5.00000E+02 -65942+25 . 105*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,NM;0) 5.00000E+02 -65348+25 . 105*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,SM;0) 5.00000E+02 -67856+25 . 105*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,YM;0) 5.00000E+02 -67419+25 . 105*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,ZM;0) 5.00000E+02 -87987+29 . 496*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,BM,NM;0) 5.00000E+02 +21425-2 . 5*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,BM,NM;1) 5.00000E+02 -19600+8*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,BM,SM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,BM,YM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,BM,ZM;0) 5.00000E+02 - 12060+11 . 16*T ; 6.00000E+03 N 

REF : 0 ! 

PARAMETER G (C , BM, YM, ZM) 300 +70000; 6000 N REF: 0 ! 

PARAMETER G(C,BM,ZM;1) 5.00000E+02 +13800+5 . 38*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,NM,SM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,NM,YM;0) 5.00000E+02 +21425-2 . 5*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G (C,NM, YM; 1) 5.00000E+02 +19600-8*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,NM,ZM;0) 5.00000E+02 -52400+27+T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,NM,ZM;1) 5.00000E+02 +24400+1 . 5*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,SM,YM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,SM,ZM;0) 5.00000E+02 -12000+11*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,SM,ZM;1) 5.00000E+02 +14500+6*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(C,YM,ZM;0) 5.00000E+02 - 12060+11 . 156*T; 6.00000E+03 N 

REF : 0 ! 

PARAMETER G(C,YM,ZM;1) 5.00000E+02 +13784+5 . 379*T; 6.00000E+03 N 

REF : 0 ! 


PHASE DELTA % 2 .571 .429 i 

CONSTITUENT DELTA : BM , YM : ZM : ! 

PARAMETER G (DELTA, YM : ZM; 0) 5.00000E+02 -79536+24 . 72*T; 6.00000E+03 N 

REF : 0 !. 

PARAMETER G (DELTA, BM: ZM; 0) 300 -79350+25*T; 6000 N REF : 0 ! 

PARAMETER G (DELTA, BM, YM : ZM ; 0 ) 300 +10200; 6000 N REF: 0 ! 
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PHASE H % 1 1.0 ! 

CONSTITUENT H : BM, NM, SM, YM, ZM : ! 

PARAMETER G(H,BM;0) 2.98150E+02 -55919+20 . 92*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,NM;0) 2.98150E + 02 - 75697 + 29 . 29*T ; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,SM;0) 2.98150E+02 -54000+30*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,YM;0) 5.00000E+02 - 5673 5+2 0 . 92 *T ; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,ZM; 0) 5.00000E+02 - 53 973 + 25 . 104 *T ; 6.00000E+03 N REF : 0 ! 
PARAMETER G (H, BM, NM; 0 ) 5.00000E+02 +12750+1 . 5*T; 6.00000E+03 N REF : 0 ! 
PARAMETER G { H , BM , NM ; 1 ) 5.00000E+02 -1000-T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,BM,SM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G ( H , BM , YM ; 0 ) 5.00000E+02 +20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,BM,ZM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,NM,SM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G ( H , NM , YM ; 0 ) 5.00000E+02 +11265+1 . 5*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G ( H , NM , YM ; 1 ) 5.00000E+02 +1000+T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,NM,ZM;0) 5.00000E+02 35000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,SM,YM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,SM,ZM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(H,YM,ZM;0) 5.00000E+02 50420; 6.00000E+03 N REF : 0 i 


PHASE LIQ % 1 1.0 ! 

CONSTITUENT LIQ : BM, NM, SM, YM, ZM : ! 

PARA G (LIQ, BM; 0) 298.15 0; 6000 N! 

PARA G (LIQ , NM ; 0 ) 298.15 0; 6000 N! 

PARA G (LIQ, SM;0) 298.15 0; 6000 N! 

PARA G {LIQ, YM; 0) 298.15 0; 6000 N! 

PARA G (LIQ , ZM ; 0 ) 298.15 0; 6000 N! 

PARAMETER G (LIQ ; BM , NM ; 0 ) 5.00000E+02 -200; 6.00000E+03 N REF : 0 ! 

PARAMETER G (LIQ , BM, SM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (LIQ, BM, YM; 0) 5.00000E+02 -10000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (LIQ, BM, ZM; 0 ) 5.00000E+02 -183750+72 . 4*T; 6.00000E+03 N 

REF : 0 ! 

PARAMETER G (LIQ, BM,ZM;1) 5.0000QE+02 +48700 - 9 . 48*T; 6.00000E+03 N 

REF : 0 ! 

PARAMETER G (LIQ, NM, SM; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (LIQ , NM, YM ; 0 ) 5.00000E+02 -200; 6.00000E+03 N REF : 0 ! 

PARAMETER G ( LIQ , NM , ZM ; 0 ) 5.00000E+02 -20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (LIQ , NM , ZM ; 1 ) 500 -15000; 6000 N REF: 0 ! 

PARAMETER G (LIQ , SM , YM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (LIQ, SM,ZM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 I 

PARAMETER G (LIQ, YM, ZM; 0) 5.00000E+02 - 183751+72 . 4 *T ; 6.00000E+03 N 

REF : 0 ! 

PARAMETER G (LIQ, YM, ZM; 1} 5.00000E + 02 +48 733 - 9 . 476*T ; 6.00000E+03 N 

REF : 0 l 


PHASE M % 1 1.0 ! 

CONSTITUENT M : BM, NM, SM, YM, ZM : ! 

PARAMETER G(M,BM;0) 5.00000E+02 -14000+29 . 915*T; 6.00000E+03 N REF : 0 

PARAMETER G(M,NM;0) 5.00000E+02 -14000+29 . 915*T; 6.00000E+03 N REF : 0 

PARAMETER G(M,SM;0) 5.00000E+02 -14000+29 . 915*T; 6.00000E+03 N REF : 0 

PARAMETER G(M,YM;0) 5.00000E+02 -14698+29 . 915*T; 6.00000E+03 N REF : 0 

PARAMETER G(M,ZM;0) 5.00000E+02 - 99979+3 5 . 9*T; 6.00000E+03 N REF : 0 ! 
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PARAMETER G(M,BM,NM;0) 5.00000E+02 -5000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (M, BM, NM; 1) 5.00000E+02 100; 6.00000E+03 N REF : 0 ! 
PARAMETER G(M,BM,SM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 

PARAMETER G ( M , BM , YM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 

PARAMETER G(M,BM,ZM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 

PARAMETER G(M,NM,SM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 

PARAMETER G ( M , NM , YM ; 0 ) 5.00000E+02 -50000; 6.00000E+03 N REF : 0 
PARAMETER G (M, NM, YM; 1) 5.00000E+02 1000; 6.00000E+03 N REF : 0 ! 
PARAMETER G(M,NM,ZM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 

PARAMETER G(M,SM,YM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 

PARAMETER G(M,SM,ZM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 

PARAMETER G(M,YM,ZM;0) 5.00000E+02 -58223+98 . 126*T; 6.00000E+03 

REF : 0 ! 


i 


N 


PHASE MN % 1 1.0 ! 

CONSTITUENT MN : BM , NM , SM , YM , ZM : 


PARAMETER G ( MN , BM ; 0 ) 5.00000E+02 -45000+30*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G ( MN , NM ; 0 ) 5.00000E+02 -60000+30*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(MN,SM;0) 5.00000E+02 -60000+30*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G ( MN , YM ; 0 ) 5.00000E + 02 -45000 + 30*T; 6.00000E + 03 N REF : 0 ! 

PARAMETER G(MN,ZM;0) 5.00000E+02 -45000+30*T; 6.00000E+03 N REF : 0 J 

PARAMETER G (MN , BM , NM ; 0 ) 5.00000E+02 -160000+48*T; 6.00000E+03 N REF : 0 

PARAMETER G (MN , BM , NM ; 1 ) 5.00000E+02 +120000-15*T; 6.00000E+03 N REF : 0 

PARAMETER G (MN , BM , SM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (MN , BM , YM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (MN , BM , ZM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (MN , NM , SM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (MN, NM, YM; 0 ) 5.00000E+02 - 160000+50*T ; 6.00000E+03 N REF : 0 

PARAMETER G (MN , NM , YM ; 1 ) 5.00000E+02 - 120000+15*T ; 6.00000E+03 N REF : 0 

PARAMETER G (MN , NM , ZM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (MN , SM , YM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (MN , SM , ZM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (MN , YM , ZM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 


; 

i 


i 

i 


PHASE N % 1 1.0 ! 

CONSTITUENT N : BM , NM , SM , YM , ZM : ! 

PARAMETER G(N,BM;0) 5.00000E+02 -50000+30*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(N,NM;0) 5.00000E+02 - 78 00 0 + 3 0 . 29*T ; 6.00000E+03 N REF : 0 J 

PARAMETER G(N,SM;0) 5.00000E+02 -50000+30*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(N,YM;0) 5.00000E+02 -50000+30*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G(N,ZM;0) 5.00000E+02 -50000+30*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G (N, BM, NM; 0 ) 5.00000E+02 -15200-7*T; 6000 N REF : 0 ! 

PARAMETER G (N, BM, NM; 1) 500 - 13 000 - 6 . 7*T ; 6000 N REF : 0 ! 

PARAMETER G(N,BM,SM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G ( N , BM , YM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(N,BM,ZM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 J 

PARAMETER G(N,NM,SM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (N, NM, YM; 0 ) 500 -15000-7*T; 6.00000E+03 N REF : 0 ! 

PARAMETER G (N, NM, YM; 1) 500 13000+5*T; 6000 N REF: 0 ! 

PARAMETER G (N , NM , ZM ; 0 ) 5.00000E+02 50000; 6.00000E+03 N REF : 0 ! 

PARAMETER G (N , SM , YM ; 0 ) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 

PARAMETER G(N,SM,ZM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 i 

PARAMETER G(N,YM,ZM;0) 5.00000E+02 20000; 6.00000E+03 N REF : 0 ! 
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PHASE P % 2 1 1 i 

CONSTITUENT P : NM : ZM : l 

PARAMETER G(P,NM:ZM;0) 5. OOOOOE+02 - 180450+57*T ; 6. OOOOOE+03 NREF:0 I 


PHASE RB % 2 4 31 

CONSTITUENT RB : SM : ZM : J 

PARAMETER G (RB , SM : ZM; 0 ) 5. 00000E+02 - 5572 0 0+171 . 5*T ; 6. OOOOOE+03 N 

REF : 0 l 


PHASE RH % 2 2 51 

CONSTITUENT RH : SM : ZM : I 

PARAMETER G (RH, SM : ZM ; 0 ) 5. 00000E+02 - 620000+185*T ; 6.00000E+03 N 

REF : 0 i 


PHASE RM % 2 2 71 

CONSTITUENT RM : SM : ZM : l 

PARAMETER G (RM , SM : ZM ; 0 ) 5.00000E+02 -822800+258*T; 6.00000E+03 N 

REF : 0 l 


PHASE T % 1 1.0 ! 

CONSTITUENT T : BM , NM , SM , YM , ZM : ! 


PARAMETER 

G ( T , BM ; 0 ) 

5 . 00000E+02 

-35000+26. 462*T 

6 . OOOOOE+03 

N 

REF : 0 J 

PARAMETER 

G ( T , NM ; 0 ) 

5. 00000E+02 

-35000+26 . 462 *T 

6 . 00000E+03 

N 

REF : 0 

PARAMETER 

G ( T , SM ; 0 ) 

5 . OOOOOE+02 

-35000+26 .462*T 

6. OOOOOE+03 

N 

REF : 0 

PARAMETER 

G { T , YM ; 0 ) 

5. OOOOOE+02 

-35618+26 .462*T 

6. OOOOOE+03 

N 

REF : 0 l 

PARAMETER 

G { T , ZM ; 0 ) 

5. 00000E+02 

-93955+31. 755*T 

6. OOOOOE+03 

N 

REF : 0 l 


PARAMETER 

G ( T , BM , NM ; 0 ) 

5 . 00000E+02 

20000, 


6 . 00000E+03 

N 

REF : 0 J 

PARAMETER 

G { T , BM , SM ; 0 ) 

5. 00000E+02 

20000, 


6. Q0000E+03 

N 

REF : 0 ! 

PARAMETER 

G { T , BM , YM ; 0 ) 

5 . OOOOOE+02 

20000 , 


6. 00000E+03 

N 

REF : 0 J 

PARAMETER 

G (T, BM, ZM; 0) 

5 . 00000E+02 

-25800; 

6. OOOOOE+03 

N REF : 0 

PARAMETER 

G (T , NM , SM ; 0 ) 

5 . OOOOOE+02 

20000 


6 . 00000E+03 

N 

REF : 0 I 

PARAMETER 

G ( T , NM , YM ; 0 ) 

5. 00000E+02 

20000 


6. OOOOOE+03 

N 

REF : 0 l 

PARAMETER 

G (T , NM , ZM ; 0 ) 

5. 00000E + 02 

20000, 


6. 00000E+03 

N 

REF : 0 1 

PARAMETER 

G ( T , SM , YM ; 0 ) 

5 . 00000E + 02 

20000, 


6. 00000E+03 

N 

REF : 0 < 

PARAMETER 

G ( T , SM , ZM ; 0 ) 

5. OOOOOE+02 

20000, 


6. 00000E+03 

N 

REF : 0 l 

PARAMETER 

G ( T , YM , ZM ; 0 ) 

5. OOOOOE+02 

-25800; 

6. 00000E+03 

N REF : 0 


LI S T_OF_RE FEREN CE S 
NUMBER SOURCE 

i 
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APPENDIX n 


$ MACRO bmymzm-3 . tcm GENERATED ON PC/WINDOWS NT DATE 0- 3- 8 
$ Make the BM-YM-ZM ternary with a L (C, BM, YM, ZM; 0) term directly 
$ 

go da 

sw user nasaSrl 
def-sys bm ym zm 
get 

go gib 

1-p-d 

c 

ent -par 

L 

C 

bm 

ym 

zm 

0 

500 

@? Enter Value 

6000 

N 

1 -p-d 
c 

go p-3 
sp-op 

SET_MI S C I B I L I TY_GAP 
c 

2 

bm ym 

s-c n=l, t=1600, p=le5, x(ym)=.78, x(zm)=.219 

1 -c 

c-e 


N 

Y 
zm 

Y 

bm ym 

N 

N 

bm 

N 

bm 

N 
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bm 

N 

bm 

N 

bm 

N 

bm 


s-c n=l, t=1600, p=le5, x (ym) = . 78 , x(zm)=.219 
1 - c 
c-e * 

1-e 

SCREEN 

WCS 

s-a-v 1 x (ym) 

0 

1 

. 025 

s-a-v 2 x(zm) 

0 

1 

. 025 
add 
-1 

li - in-eq 
save 

bmymzm-3 

Y 

map 
post 
s-d- 1 

y 

Y 
n 

pl 

SCREEN 

ba 

read 

bmymzm- 3 . PL3 

1 -c 

s-c x (ym) =le- 3 , x(zm) = .21 

1-c 

c-e * 

c-e 

1-e 

SCREEN 

VWCS 

add 

2 

li-in-eq 

s-c x (ym) = . 05 , x(zm)=.9 

1-c 

c-e 

c-e 

1-e 

SCREEN 

VWCS 
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add 

-2 

li-in-eq 

save 

bmymzm- 3 . PL3 

y 

map 
post 
s- 1- s 
4 

s - c - t L_L 
YbOl . 5 
s-c-t 1-r 
YOl . 5 
s-c-t t 
Zr02 
s- f 
4 

0.3 

s-a- t- s 

x 

N 

X (YOl .5) 
s-a- t- s 
Y 
N 

X (Zr02) 
s- tit 

YbOl . 5 -YOl . 5 - Zr02 1600K 
s - 1 - c 
n 
Pi 

SCREEN 

dump 

bmp 

ba 

read bmymzm- 3. PL3 
li-in-eq 

s-c x(ym)=0.78, x(zm)=0.21 

1-c 

c-e * 


1-e 


add 

-1 

save 

bmymzm- 3 . PL 3 

y 

map 

post 

Pi 

SCREEN 
set- int 
ba 
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TABLE I.— AVAILABLE BINARY PHASE DIAGRAMS 



ZiA 

Y A 

YbA 

Sep, 

Ndp, 

ZrO, 


X 

X 

X 

X 

YA 



X 

X 

X 

YbA 




X 

X 

ScA 





X 

Nd.O, 







TABLE n.— ZIRCONIA AND RARE EARTH OXIDES TRANSITION TEMPERATURES 


AND HEATS OF TRANSFORMATION 


Oxide 

Phase Change 

T(trans) 

K 

A.H 

(kj/mol) 

A„S 

(J/mol-K) 

Reference 

ZrO, 

Monoclinic — » Tetragonal 

1448 

5.941 

4.103 

Kubachewski (33) 


Monoclinic — » Tetragonal 

1478 

8.075 

5.464 

Pankratz (35) 


Monoclinic — > Tetragonal 

1445 

8.4 

5.813 

IVT AN (34) 


Tetragonal — » Cubic 

2568 

5.9 

2.3 

Kubachewski (33) 


Tetragonal — > Cubic 

2620 

13.001 

4.962 

IVT AN (34) 


Cubic — > Liquid 

2953 

87 

29.492 

Kubachewski (33) 


Cubic — » Liquid 

2983 

90 

30.171 

IVT AN (34) 







YA 

Cubic — » Hexagonal 

2550 

54 

21.176 

IVT AN (34) 


Hexagonal — » Liquid 

2712 



Kubacheski (33) 


Hexagonal — » Liquid 

2712 

81 

29.867 

IVT AN (34) 







YbA 

Cubic — » Hexagonal 

2663 



IVT AN (34) 


Hexagonal — » Liquid 

2708 

130 

48.006 

IVT AN (34) 







NdA 

N phase — » Hexagonal 






Hexagonal — > Cubic 

2333 



IVT AN (34) 


Cubic — » Liquid 

2593 

125 

48.207 

IVT AN (34) 
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TABLE III — OXIDES 
AND ABBREVIATIONS 


Oxide 

Abbreviation 

ZtQ 2 

ZM 

YO,3 

YM 

YbO u 

BM 

NdO, j 

NM 

ScO J3 

SM 


TABLE IV. — MOLE FRACTION OF RE,0 3 IN RE£> r ZRO, 
AND EQUIVALENT MOLE FRACTION OF 
REO, s IN REO, ,-ZRQ, 


Mole Fraction RE 2 Oj 

Equivalent Mole Fraction REO, 5 

x(RE 3 0,) 

x(REO, 5 ) 

0.1 

0.182 

0.2 

0.333 

0.3 

0.462 

0.4 

0.571 

0.5 

0.667 

676 

0.75 

0.7 

0.824 

0.8 

0.889 

0.9 

0.947 


TABLE V.— EFFECT OFTTRST ORDER 
INTERACTION PARAMETER ON 
THE LOCATION OF THE 
MISCIBILITY GAP 


Miscibility Gap 

X 

Symmetric 

Zero 

Shifted to Right 

Positive 

Shifted to Left 

Negative 
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Figure 1 .—Two experimentally determined Y2C>3-Zr02 phase diagrams, 
(a) Stubican et al. 1981 (22, 7). (b) Pascual and Duran, 1983 (23, 9). 
Reprinted with permission of The American Ceramic Society, Post 
Office Box 6136, Westerville, OH 43086-6136. Copyright [1984] by 
the American Ceramic Society. All rights reserved. 
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Temperature, °C Temperature, 



1 800 l— ' 1 1 1 1 1 

0 20 40 60 80 100 


Zr °2 Mol % y 2°3 



Figure 2. — Experimentally determined Yb2C>3-Zr02 phase diagrams, (a) Rouanet 
(24). (b) Stubican et al. (11). Reprinted with permission of The American Ceramic 
Society, Post Office Box 6136, Westerville, OH 43086-6136. Copyright [1985] 
by the American Ceramic Society. All rights reserved. 
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Figure 3. — Experimentally determined SC2O3-Z1O2 phase diagram 
( 25 ). Reprinted with permission of The American Ceramic Society, 
Post Office Box 6136, Westerville, OH 43086-6136. Copyright 
[1993] by the American Ceramic Society. All rights reserved. 


NASA/TM— 200 1 -2 10753 


23 




Zr °2 Mol % Nd 2 0 3 

Figure 4. — Experimentally determined Nd 203 -ZK >2 phase diagram (26). 
Reprinted with permission of The American Ceramic Society, Post Office 
Box 6136, Westerville, OH 43086-6136. Copyright [1993] by the American 
Ceramic Society. All rights reserved. 
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Figure 5. — Experimentally determined Yb203-Y203 
phase diagram (27). Reprinted with permission of 
The American Ceramic Society, Post Office Box 6136, 
Westerville, OH 43086-6136. Copyright [1993] by the 
American Ceramic Society. All rights reserved. 



y 2^3 Mol % Sc 2°3 

Figure 6.— Experimentally determined SC2O3-Y2O3 
phase diagram (28). Reprinted with permission of 
The American Ceramic Society, Post Office Box 
6136, Westerville, OH 43086-6136. Copyright [1993] 
by the American Ceramic Society. All rights reservec 
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Figure 7. — Experimentally determined Nd203-Y2C>3 phase diagram 
(29). Reprinted with permission of The American Ceramic Society, 
Post Office Box 6136, Westerville, OH 43086-6136. Copyright 
[1993] by the American Ceramic Society. AH rights reserved. 



Yfc>20 3 Mol % Sc 2 0 3 


Figure 8.— Experimentally determined Sc203-Yb203 
phase diagram (28). Reprinted with permission of 
The American Ceramic Society, Post Office Box 
6136, Westerville, OH 43086-6136. Copyright [1993] 
by the American Ceramic Society. All rights reserved. 
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Figure 9. — Experimentally determined Nd203-Yb2C>3 phase diagram (30). 
Reprinted with permission of The American Ceramic Society, Post Office 
Box 6136, Westerville, OH 43086-6136. Copyright [1993] by the American 
Ceramic Society. All rights reserved. 
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Figure 10. — Experimentally determined Nd203-Sc203 
phase diagram (28). Reprinted with permission of The 
American Ceramic Society, Post Office Box 61 36, 
Westerville, OH 43086-6136. Copyright [1993] by the 
American Ceramic Society. All rights reserved. 


NASA/TM — 2001-210753 


27 



^ss\ 1 200 °C 


Fss + Zr 3 (Y, Yb) 4 0-|2 


I I \ 

Zr3(Y, Yb) 4 0| 2 + Css 

/ I \ 


// I '\ 
Tbs + Fss 

! ! \ 


1400 °C 


Fss + Zr 3 {Y > Yb) 4 0 12 


iS^^SS 

+ Zr 3 (Y,Yb)40 1 2 


\ \ ZT3(Y,Yb)40|2 \ v \ 

\ C^ + Zr 3 (Y,Yb) 4 0 12 \ 

\ \ \ \ \ 


Figure 1 1 .—Experimental determined Y203-Yb203-Zr02 phase diagram (31) at (a) 1 200 °C, (b) 1 400 °C, 
and (c) 1650 °C. Reprinted with permission of The American Ceramic Society, Post Office Box 6136, 
Westerville, OH 43086-6136. Copyright [1993] by the American Ceramic Society. All rights reserved. 
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Figure 1 1 . — Concluded, (c) 1650 °C. Reprinted with permission of The American 
Ceramic Society, Post Office Box 6136, Westerville, OH 43086-6136. Copyright 
[1993] by the American Ceramic Society. All rights reserved. 



xi 

Figure 12.— Values of G con fj g + G eX cess versus T illustrating the effect of the 
interaction parameters on the symmetry and location of the miscibility gap. 
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Zr °2 xfYO!,^ Y0 1.5 


Figure 1 3. — Calculated YOj g-ZrOg phase diagram. Compare to Figure 1 (a). 
The starting point for the calculation is x(YM) = 0.5, T = 3500 K. 
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x<ybOi' 5 ) YbOi. 


Figure 14.— Calculated Yb0i.5-Zr02 phase diagram. Compare to Figure 2. 
The starting point for the calculation is x(ZM) = 0.5, T = 3500 K, followed 
x(ZM) = 0.2, T = 2000 K. 


N AS A7TM— 200 1 -2 1 0753 


31 




Zr0 2 xtScO-,. 5 ) ScO-| 

Figure 15. — Calculated SCO 1 . 5 -Z 1 O 2 phase diagram. Compare to Figure 3. 
The starting point for the calculation is x(ZM) = 0.5, T = 3500 K. 
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Zr0 2 xfNdOt.s) Nd0 1-5 

Figure 16. — Calculated Nd0i.5*Zr02 phase diagram. Compare to Figure 4. 

The starting point for the calculation is x(ZM) = 0.1 , T = 1800 K, initially 
mapping up with temperature, followed by x(ZM) = 0.99, T = 1000 K, 
initially mapping up with temperature, followed by x(ZM) = 0.99, T = 3100 K, 
initially mapping up with temperature, and finally x(ZM) = 0.1, T = 2300 K. 
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Figure 17. — Calculated YbO-i.s-YO-j.s phase diagram. Compare to Figure 
5. The starting point for this calculation is x(YM) = 0.001 , T = 2500 K, 
map with temperature increasing first. 


NASA/TM — 200 1 -2 1 0753 


34 



Temperat 


Liquid 


l 
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ScO-j .5 x(Y0 1-5 ) YOv 


Figure 1 8.— Calculated SCO 1 . 5 -YO 1.5 phase diagram. Compare to Figure 
6 . The starting point for this calculation is x(SM) = 0.001 , T = 2500 K, 
map with temperature increasing first. 
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Temperature, 



Nd0 1-5 xfYCH.s) YOu 

Figure 1 9.— Calculated NdO-i.s-YOi.s phase diagram. Compare to Figure 
7. The starting point for this calculation is x(YM) = 0.01 , T = 3000 K, 
followed by x(YM) = 0.99, T = 1 500 K. 
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Temperature, 



0.0 0.2 0.4 0.6 0.8 1.0 

YCH.s x (ScO-1.5) ScO 


Figure 20.— Calculated ScO^-YO^ phase diagram. Compare to Figure 
8. The starting point for this calculation is x(SM) = 0.001 , T = 2500 K, 
map with temperature increasing first. 


NASA/TM— 2001-210753 


37 



Temperature, 



NtiOi .5 x(YbOi.s) YbOi.j 

Figure 21 -Calculated NdOi.s-YbO-j.s phase diagram. Compare to Figure 
9. The smarting point for this calculation is x(BM) = 0.99, T = 3000 K, 
followe by x(BM) = 0.1, T = 1500 K. 
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Temperature, 



NdO -|.5 x (ScO-| . 5 ) ScO, ,j 

Figure 22. — Calculated ScOi . 5 -NdOi .5 phase diagram. Compare to Figure 
10. The starting point for this calculation is x(SM) = 0.001 , T = 3000 K, 
followed by x{SM) = 0.001, T = 1000 K. 


NASA/TM— 2001-2 10753 


39 
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x (Y0 1<5 ) 
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Figure 23. — Calculated activities across the YOi.s-ZrC^ 
phase diagram at 2773 K (line) and experimental data 
(points) of Belov and Semenov et al. (39). This line was 
generated using a starting point of x(ZM) = 0.8 and 
then a step with options command, (a) a (YO 1 . 5 ). 

(b) a (ZrOa). 
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x (ScOi .5) 



0.0 0.2 0.4 0.6 0.8 1.0 

x (ScO-i.s) 

Figure 24.— Calculated activities across the ScO*|.5-Zr02 
phase diagram at 2600 K (line) and experimental data 
(points) of Belov and Semenov et al. (40). This line was 
generated using a starting point of x(ZM) = 0.8 and 
then a step with options command, (a) a (ScOi.5). 

(b) a (Zr02). 
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Tet ZrC>2 + cub ZrC>2 


Zr0 2 



YbOi 5 0 0.2 0.4 0.6 0.8 1.0 yo 1-5 

x (Y0 1-5 ) 

Zr0 2 



YbOi.s 0 0.2 0.4 0.6 0.8 1.0 YO^ 

x (YO1.5) 

Figure 25. — Calculated ternary section for the YbOi 5-YO1 5 -Zr0 2 
system at 1600 K. (a) L(C, BM, YM, ZM; 0) = +10 000 insufficient 
to form a continuous two phase region, (b) L(C, BM, YM, ZM; 0) = 
+100 000 continuous two phase region. 
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Temperature, °C 


Figure 27. — Calculated phase fraction diagram for 
0.05YbQ 1 .5-0.04YO-1 .5-O.OI ScO, . 5 -0.9ZrO 2 . 



Figure 28. — Calculated isopleth for 
Yb0 1 ,5-Y0 1 .s-NdOi . 5 -Zr0 2 with 
x(Yb0i.5) = x(Sc0 1 .5) = 0.02. 
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